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Abstract 

We consider real eigen-functions of the Schrodinger operator in 2-d. The nodal lines of separable 
systems form a regular grid, and the number of nodal crossings equals the number of nodal domains. 
In contrast, for wave functions of non integrable systems nodal intersections are rare, and for 
random waves, the expected number of intersections in any finite area vanishes. However, nodal 
lines display characteristic avoided crossings which we study in the present work. We define a 
measure for the avoidance range and compute its distribution for the random waves ensemble. We 
show that the avoidance range distribution of wave functions of chaotic systems follow the expected 
random wave distributions, whereas for wave functions of classically integrable but quantum non- 
separable wave functions, the distribution is quite different. Thus, the study of the avoidance 
distribution provides more support to the conjecture that nodal structures of chaotic systems are 
reproduced by the predictions of the random waves ensemble. 
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I. INTRODUCTION 



The morphology of the nodal sets of wave functions depends crucially on whether the 
underlying classical dynamics is integrable or chaotic. This was first proposed in [l|] and 
was followed by the study of various features of the nodal lines, such as the distribution 
of its curvatures [0]. Recently, the counting statistics of nodal domains for integrable and 
chaotic systems were investigated in ||, and it was shown that in the chaotic case, the 
statistics follow the predictions derived by assuming that the wave functions are random 
superpositions of plane waves [BL Local effects due to boundary conditions, and the 
corresponding modifications of the random wave ensembles were also studied R [7J . 

The interest in the properties of the nodal set is not confined to the physics literature 
only. Most of the mathematics literature on this subject is concerned with solutions of the 
Helmholtz equation in the interior of compact domains in R 2 with Dirichlet or Neumann 
boundary conditions. Courant ||, |] and later Pleijel [JlO|] pioneered these studies, and 
computed an upper bound on the number of nodal domains. Krahn Jll| provided a lower 
bound for the area of a nodal domain. Other authors gave estimates on e.g., the length of 
the nodal sets fl2|| , and addressed the general properties of the nodal network [[13]. 



Uhlenbeck's theorem [|14|] states that the nodal lines of "generic" wave functions do not 
intersect (this statement will be further discussed in the next section). An important class 
of exceptions of this rule, are the eigen-functions of separable systems, where the nodal 
lines form a grid, and the number of intersections equals the number of nodal domains. 
The purpose of the present work is to provide a quantitative measure of the degree by which 
nodal lines avoid each other. We shall associate the avoidance range to each avoided crossing, 
and will compute their distribution for the random waves ensemble. The avoidance range 
vanishes for an intersection, and therefore, the avoidance range distribution for separable 
functions is proportional to a Dirac 5 at zero. We also show that eigen-functions of classically 
chaotic systems display an avoidance range distribution which is similar to the one obtained 
for a random waves ensemble, while for an intermediate system, the avoidance distribution 
differs substantially from either of the extreme distributions mentioned above. The nodal 
structure for a few representative wave functions can be seen in figure ([3]). 

The wave function of a rectangle (a) displays a perpendicular grid of nodal lines typical 
of separable systems. An equilateral triangle (b) is classically integrable, but not separable, 
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FIG. 1: Eigen- functions of 2-d billiards: (a) rectangle (separable dynamics); (b) equilateral triangle 
(integrable, non-separable dynamics); (c) Sinai (chaotic dynamics). 

and the wave function shows a few crossings and avoided crossings of the nodal lines. The 
domain in (c) is a Sinai billiard where the classical dynamics is chaotic. With the exception 
of the boundary crossings, the nodal set displays a characteristic set of avoided crossings. 

The paper is organized in the following way: We shall begin the next section by discussing 
the nodal structure of general solutions of the Helmholtz equation. We shall then identify 
avoided crossings of the nodal lines, and define the corresponding avoidance range. The 
avoidance range distribution will be written down explicitly. In section ( |III| ) we shall define 
the random wave ensemble and compute explicitly the expected distribution of the avoidance 
ranges. The resulting expression will be compared with the avoidance distributions obtained 
numerically for high lying eigen-functions of a chaotic domain. For Dirichlet problems, the 
boundary belongs to the nodal set, where nodal intersections occur at a density which is 
approximately 2 intersections per wavelength. This property cannot be reproduced by the 
uniform random wave ensembles, which are expected to account for bulk properties. This 
problem, in the context of the avoided crossing distribution, is discussed at the end of section 

Q 



3 



II. AVOIDED CROSSINGS OF NODAL LINES AND THE AVOIDANCE RANGE 



We consider real solutions of the Helmholtz equation in a domain Q e M 2 which can be 
expressed as linear superpositions of regular solutions of the equation in the entire plane. 
They can be expressed in terms of either plane or cylindrical waves, and the two repre- 
sentations are equivalent. In the plane waves representation the wave function is written 

as 

oo 

tt(r) = a ™ ^ W 

n=— oo 

with wave vectors (k n = — k_ n , |k n | 2 = k 2 ) directed at an angle 8 n and coefficients which 
ensure that the series (|l|) converges absolutely (we require a* = a_„ to render \l/(r) real). 
Using the expansion of a plane wave in cylindrical coordinates we get 

(r) = a, Ji{kr) e ue , (2) 

l=— oo 

with 



a>i 



CL n e iWn , and a ; * = (-l) l a_i . (3) 



The infinite sums in both (JJ) and (0) can be truncated, when one considers domains of 
finite area. Semiclassical arguments show that the necessary number of terms is L « Ck/n, 
where £ is the perimeter of the domain. For reasons which will become clear in the sequel, 
we prefer to use the cylindrical wave representation in the present work. 

The expansion ([2]) refers to a particular choice of the origin. Using Graf's addition 
theorem |T^] the origin can be shifted to r, and the wave function retains its form, 

v&(r + p) = «m(r) J m (kp) e im * , (4) 

m 

where 



a, 



(r) = £>,((>) J^Ukr) e lW , (5) 
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with ai(0) = ai, and the angle is measured from the direction defined by r. The translated 
coefficients a m (r) = (3 m {r) + ry m (r) are related to the wave-function and its derivatives 
computed at the point r: 




+ |ai(r)| ( ) cos(0 + 1 ) + ^|a 2 (r)| 1^1 cos 2(0 + 2 ) 



A)(r) = *(r) (6) 

A(r) = ±* P (r) ; 7 i(r) = -^(r) 

2 2 
/3 2 (r) = ^(r) + — ^ rr (r) ; 72 (r) = — (* e (r) - r^(r)) 

<Mr) = 4rft(r)~ (&(r) + /3 (r)) 

In the close vicinity of r, where kp < 1, and to second order in /cp 



*(r + p) « /5o(r)(l-(T) ) (7) 

2 

2 / 

where 0; are the phases of ai(r). If /3q(t) = 0, r is a nodal point. It cannot be an isolated 
zero since the second term vanishes on the line segment through r which is oriented at the 
direction | — 0i- Hence, the nodal set consists of lines. Two nodal lines intersect at r if 
both Po(r) = and ai(r) = 0, while a 2 (r) ^ 0. The intersection is perpendicular since 
cos2(0 + 2 ) vanishes along two perpendicular lines which intersect at r. For the time being 
we shall continue the discussion assuming that a 2 (r) ^ 0. The more general case will be 
commented on at the end of this section. 

An avoided crossing occurs at r when a!i(r) = and |a 2 (r)| > |/3o(r)| > 0. In other 
words, when r is a saddle point of the wave function. This can be easily seen by writing the 
equation of the zero set of (|7|) in terms of the local coordinates p = (£,T]), 

, = P 2 / A)(r)-Mr)| \ 2 / /? (r) + |q 2 (r)| \ 

4 V A)(r) ) +V V &(r) ) ' U 
This is a hyperbola (ellipse) if |a 2 (r)| is larger (smaller) than |/3 (r)|. At an avoided crossing, 

the scaled distance between the two branches is 
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This is the avoidance range associated with the avoided crossing at r. 
A few comments are in order: 

i. At a nodal intersection z = 0. 

ii. At a saddle point | ar 2 (r) | > |/3 (r)|, hence z < 2\/2. 

Hi. An equivalent expression for z(r) in terms of \I/(r) and its Cartesian derivatives reads 



z(r) = 4 



k 2 m 



(10) 



^01*1 + ^4*%, + (*„-* w y 

iv. Consider an elliptic critical point of the wave function. In the quadratic approximation 
the area of the elliptic nodal domain is 



A = 47rAr 



ft 



P 2 -\a 2 \ 2 



(11) 



This area is always larger than 4irk~ 2 . Krahn's theorem |TTJ gives j'q 1 7r/c~ 2 as the maximal 
lower bound to the area of any nodal domain, where jo,i ~ 2.405 is the first zero of the 
Bessel function Jq(x). The lower bound 47rfc -2 is smaller but not very far from Krahn's 
exact value and thus consistent. 

So far we considered the intersections of two nodal lines. However, higher order intersec- 
tions may occur. In general, if the first non vanishing coefficient at r is a q , then r is a nodal 
point of order q, where q nodal lines intersect at angles -. The higher q, the more rare are 
the intersections, since more conditions are to be satisfied by the coefficients. This explains 



Uhlenbeck's theorem [H] mention above. From now on we shall discuss the most com- 
mon intersections with q — 2, and comment about the higher order intersections whenever 
necessary. 

Up to now we discussed individual avoided crossings, and defined the associated avoidance 
ranges. In the next section we shall consider the distribution of the avoidance ranges of a 
wave function in the domain of its definition, and compute its mean for random waves 
ensembles. 
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III. AVOIDANCE RANGE DISTRIBUTIONS 



The number of the critical points of \l/(r) in the domain Q is given by 

N c = jjdr d9 5(*r(r)) S (~**(r)) \J(r)\ , (12) 
with J(r) = i(^ rr (r)^ ee (r) - ^(r)) the Jacobian. Using (0) we get 

N c = ^ I rdr d9 *(ft(r)) 5( 7 i(r)) ||a 2 (r)| 2 - /3 2 (r)| . (13) 
To count the saddle points, we add the restriction |a 2 (r)| 2 — /5o( r ) > and obtain 

iV s = ^J u rdr d6 *(ft(r)) 5( 7l (r)) (|a 2 (r)| 2 - /3 2 (r)) 6(|« 2 (r)| 2 - /3 2 (r)) . (14) 
Combining (|14[) and @, the number of avoided crossings with avoidance ranges less than z 



is 



= £ / rdrd^(/3 1 (r))5( 7l (r))(|a 2 (r)| 2 -/3 2 (r))0(|« 2 (r)| 2 -/3 2 (r))e(z-z(r)), (15) 
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and the fraction of the total number is 1{z) = T(z)/Ns- This counting function, and its 
associated density V(z) = are the distributions which characterize the avoided crossings 
of the nodal set. The range of z is [0, 2y / 2]. I(z) is normalized such that it takes the value 
1 at z = 2y/2. One can easily check that an avoided crossing of order q is counted in fllBl) 



with a multiplicity q — 1. The number (including multiplicity) of nodal crossings provides 
the value of 1(0). The effective multiplicity of nodal crossings of the boundary is reduced 
by a factor one half. 

Numerically computed T{z) are shown in figure [| for the three wave functions plotted in 
figure 0. The computed 1{z) do not take into account the nodal crossings of the boundary. 
For the separable billiard all the saddle points have zero avoidance range, and the function 
Z(z) is trivially equal to one in all the z range. In the triangular billiard there are still 
some saddle points with zero avoidance range (nodal crossings), but most of them have a 
finite avoidance range. For the chaotic billiard, all the calculated saddle points have a finite 
avoidance range. 
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FIG. 2: Normalized cumulative histograms of the avoidance ranges I(z) for the three wave functions 
plotted in figure 1: rectangle (dash-point line); triangle (dashed line); Sinai (full line). 

A. Avoidance range distribution for the random waves ensemble 

One of the main goals of the present work is to show that the properties of the nodal 
set of chaotic billiards, as detected by the distribution of avoidance ranges, are reproduced 
by the distributions computed for the random wave ensembles. Because wave functions are 
subject to boundary conditions, it is expected that the predictions of the isotropic random 
wave ensemble used e.g. in [H ||, are relevant only to the bulk of the domain, and will do 
poorly in the A = (^) vicinity of the boundary. This was observed and discussed in |J. 

The isotropic random wave ensemble is the ensemble of wave functions (0) where the real 
parameters f3i and 7/ are independent, identically distributed random Gaussian variables 
with zero mean and unit variance for all |/| > 1. Because 70 = the variance of f3o is 
twice that of all the others. The ensemble average of a function / will be denoted by (/). 
The local coefficients cty(r) = $(r) -M7z(r) were derived from the original ones by a unitary 
transformation. Hence, they are also independent identically distributed Gaussian variables, 
and the ensemble averages of the number of critical points (|12D , the number of saddle points 

8 



(|14]) and the mean distribution of avoidance ranges (|X5|) can be computed by considering 
Gaussian integrations with respect to the variables @q, fli, 71, and 72. A straightforward 
integration gives, 

P IOI A- 2 IOI 

where is the area of the domain. The average number of saddle points is a half of the 
number of critical points. The other half are the points were the wave function has either a 
minimum or a maximum. 

The mean number of saddle points with an avoidance range less than or equal to z can 
be also computed, 

,~, „ k 2 \Q\ 3 z 2 (16 -z 2 ) 2 r , . 

The normalization gives trivially the mean number of saddle points in the ensemble. There- 
fore we can define the probability 

X (z) {i{z)) - 3 v^ 2 ( 16 -^) 2 <z<2v/2 (18) 

T ^- (Z) - W" (512-64^ + 3^)3 /2 ' ° <Z<2V2 ' (18) 

with a corresponding density 

_ , , 6144^ z (8 - 2 2 ) (16 - z 2 ) „ n y- 

These results served to test the conjecture that the bulk properties of the nodal sets 
of chaotic wave functions are reproduced by the predictions of the isotropic random wave 
ensemble. For this purpose we computed numerically the first 2400 eigen-functions of the 
Sinai billiard shown in figure p].c. For each wave-function the critical points were found by 
a numerical search. The saddles which correspond to boundary intersections were excluded. 
The avoidance range ( |10|) was computed for each saddle. In figure |3| we plot the cumulative 



histogram of the avoidance ranges for the 2000th eigen-state in full line (1102 saddle points), 
compared with formula ( [P7D in dashed line. Even for a single eigen-state the agreement is 
very good. Averaging the avoidance range distributions over a group of neighboring eigen- 
states the numerical histogram and the theoretical curve approximately coincide. The inset 
in figure ^| shows the differences between the random wave prediction and the mean avoidance 
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0.002 




FIG. 3: Counting function X{z) for the n = 2000 eigen-state of the Sinai type billiard (full thin 
line), compared with the random wave prediction (|l8|) (dash line). Inset: difference between the 
average counting function over the n = (1800, 2000) eigen-states and the random wave prediction 
(dark line); the same for n = (2200, 2400) eigen-states (light line). 

range distributions computed for two groups of eigen-states. The difference is small, and 
shows no systematic deviations. 

Another distribution which we compared to the prediction of the random waves ensemble 
is the number of saddle points. Normalizing the number of saddle points by the predic- 
tion of the isotropic random wave ensemble, we observe as k increases that the numerical 
computation approach the predicted value (see figure [|). 

The systematic deviation observed at finite k is due to boundary effects. Following || 
we compute the effect of an infinite straight Dirichlet line, which is reproduced by the wave 
ensemble 

oo 

tf(r) =2^c„ sm(n6) J n (kr) , (20) 

n=l 

where the c n are real coefficients taken as independent random Gaussian distributed vari- 
ables. 
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Not entering into the details of the calculation, the computed mean density of saddle 
points approach the bulk expression ( |T6| ) as the distance from the Dirichlet line increases. 
The integrated density in the perpendicular direction shows a global deficiency of saddle 
points relative to the bulk value. It diverges logarithmically as a function of the distance 
from the Dirichlet line 0. 

We compared the results of the ensemble fl2"UD to the number of saddle points counted 
for highly exited eigen-states of a chaotic billiard with sufficiently smooth boundaries. The 
integral along the boundary multiplies the density by the perimeter of the billiard C The 
perpendicular integral must be truncated in view of the logarithmic divergence mentioned 
above. A sensible choice of the truncation distance is R = ^/Jfi[/2. The resulting estimate 
for the mean number of saddle points is 

^ k*\{l\ - kCja^ogjkR) + a 2 ) 
[ 5/ ~ 4vrx/3 ' 1 ] 

with o\ « 0.014 and o~i ~ 2.0. The deficiency of saddle points is explained by the effect of 

the Dirichlet boundaries that affects the statistics. The dashed line in figure [| represents the 

equation ( pTD and the agreement is definitely improving. However, the domain of low values 

k < 500 is not in complete agreement. This deficiency can be associated with the corners 

of the billiard, to its curvature, or to the finiteness of the Dirichlet line that the ensemble 

(p0|) can not reproduce. These possible causes will be studied elsewhere. 

In conclusion we can say that the properties of the nodal set of chaotic wave functions 

which were investigated in this work, are very well reproduced by the isotropic random waves 

ensemble in the semiclassical domain. These findings are consistent with previous works on 

the subject, and add support to the random waves conjecture. 
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FIG. 4: Number of saddle points for the first 2400 eigen-states of the Sinai billiard, normalized 
by the isotropic random wave prediction k 2 /(4iry3), as a function of the wave number k. Dashed 
line: equation (f2~I|) for the anisotropic random wave ensemble (|20|). 
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